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Iron Line Emission from a High Temperature 
Plasma in Cas A 
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R. E. ROTHSCHILD, P. J. SERLEMITSOS, J. H. SWANK^ 

NASA/ Goddard Space Flight Center 
Greenbelt, Maryland 20771 

ABSTRACT 

The x-ray spectrum of Cassiopoia A was observed for several 

days by the GSFC proportional counter experiment on board OSO-8. 

The high energy (> 5 keV) data are fit well by a thermal spectrum 
+ 9 

with kT “ 3.9 _ ^ keV. A narrow iron line which is predicted by the 

thermal model is also observed, centered at 6.66 keV with an 

equivalent width of 1270 eV. The low energy (2-5 keV) data show 

an excess over the high temperature component which is consistent with 

the presence of an additional low temperature thermal component with 

22 -2 

kT ^ 0.7 keV and Njj 10 atoms cm . Iron abundance in the source 
relative to normal cosmic abundance is discussed, as in the relation 
of this observation to shock wave and multi-component thermal models 
for supernova remnants. 


*Work supported by U. of Maryland Grant 21-002-316 
^NAS/NRC Research Associate 
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I. Introduction 

Casslopola A has the highest Intrinsic and apparent X-ray luminosity 
of all known supernova remnants except for the Crab Nebula. Optical 
(van den Bergh and Dodd 1970), radio (Rosenberg 1970 a, b) , and X-ray 
(Fabian 1973) observations have shown this most recent supernova 

(100-400 years old) to be an extended source. In this Letter we report 
on the measurement of the X-ray spectrum of Cas A with a proportional 
counter on board Orbiting Solar Observatory 8 (OSO-8) . 

II. Experiment 

OSO-8 was launched June 21, 1975. It carries seven experiments. 

The two prime solar experiments are mounted on the spacecraft sail 
and five others on the wheel, the spinning portion of the spacecraft. 

One of this latter group is the GSFC Cosmic X-ray Spectroscopy experi- 
ment (CXS) . It consists of three multi-anode, multi-layer gas pro- 
portional counters in or near alignment with the spacecraft spin axis. 

The detector relevant to this observation is a xenon-filled 
counter featuring an additional gas volume of neon-propane adjacent to 
the collimator to guard against precipitating electrons entering via 
the collimator. This detector has a 5 degree circular field of view 
centered about the positive spin axis. The view is periodically 
occulted by the pointed solar instruments and by two occultation 
shields mounted under the spacecraft sail. These are designed to pro- 
vide a measure of the intrinsic detector background. Data from this 
detector are binned by onboard logic that distinguishes whether the 
field of view is open or occulted. Other relevant features include a 
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net area of 237 cm » a transmission window of two .001" mylar sheets 
in addition to the neon-propane gas in the anticoincidence volume, and 
an onboard Am^^^ calibration source. 

The observation was carried out during the period July 17-26. 

The results presented here are from "Quick Look" data corresponding 
to about 10% of all data. Background measurements were made just 
before as well as after Cas A was in the field of view. The background 
in this region of the sky was found to be very stable. 

III. Results 

The excellent statistics allow a multi-parameter fit to the Cas A 

spectrum from which meaningful model-building can proceed. First the 

high energy (> 7 keV) points alone were fit, agreeing equally well with 

a power law with number index 4.5 +.5 and a thermal spectrum with 

kT = 3.2 keV, where the errors quoted refer to acceptability at 

90% confidence. These confidence intervals are calculated by allowing 

2 

the minimum acceptable value for X to vary by an amount consistent 
with the number of parameters in the model (Lampton ^ 1975) . The 

analytical form used to represent a thermal spectrum is 

HN 

^ = C g(E,kT) exp (-E/kT)/E (1) 

-0.4 

with C a normalization constant and g(E,kT) * (E/kT) * an approximate 
form for the Gaunt factor. Tlnis spectral form (equation 1) adequately 
approximates the continuum contribution from an isothermal optically 
thin tliermal plasma calculated in various numerical models (see below) . 
The model in equation 1 can be tested against the data either by folding 
it through the detector response function (forward direction) or by 
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multiplying it by spectral dependent channel-by-channel efficiencies 

(backward direction) when good fits are obtaim:d. These two methods 

converge to the snme best fit parameters. 

Next, the energy region around 7 keV is Included in the analysis. 

An acceptable best fit is found by Including a narrow line centered at 

6.b6 keV. When the detector resolution was introduced as a free 

parameter and varied, the best fit was found at the resolution value 

previously determined by calibration (FWHM 1.1 keV at 6.66 keV) . The 

^25 ^2A 

2p -* IS transitions in thermally excited Fe and Fe would form a 

feature with these characteristics. The iron line has an equivalent width 

of 1270 eV, while the channel energy width in the CXS is ~ .6 keV. 

Accordingly, the line to continuum ratio in the channel centered at 6.7 

+ .12 

keV, where most (~ 45%) of the line photons fall, is .96 _ q^. Some 

evidence is seen for an enhancement around 8.2 keV, where lines from np -• IS 
^25 "^2A 

(n > 2) of Fe and Fe are expected at a lower level. 

The thermal spectrum and iron line, which yield an acceptable fit 

to data points above 6 keV, do not fully account for the lower energy 

data. The residuals are fit well with an additional low temperature 

22 -2 

thermal with kT ^ . 7 keV and N„ ~ 10 atoms cm . The contribution 

H 

of this component changes the best fit parameter for the high temperature 

+.9 

component, which now becomes 3.9 _ ^ keV. 

We present the data and best fit model in the first two figures. 

Figure la shows the data as points and bars representing the statistical 
errors, superimposed on a histogram of the CXS response when the input 
spectrum is the best fit model described above. A value of 23 
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wlth 19 degrees of freedom is obtained. The iron line is represented 
as a Gaussian centered at 6.66 keV with FWHM « 1.1 keV. If the FWHM 
of the Gaussian is Increased to be fm 1.1 keV the acceptable fit is lost. 
Since we were able to find convergent best fit models, accurate channel- 
by-channel efficiencies were calculated. Figure lb shows the inferred 
incident spectrum obtained by analysis in the backward direction making 
use of the efficiencies described above. 

In Figure 2 the observed fc-.ture around 6.7 keV is displayed by 
folding the best fit therma] continuum through the detector response 
function and subtracting this resultant from the data. The width and 
shape of this feature are consistent with detector broadening of a 
narrow line. Although the observed line falls in several channels, 
the analysis described above allows us to conclude that with better 
resolution, we would have observed it in only one channel (Fig. lb). 

In contrast to the complex spectrum of Cas A, Figure 3 shows the 
spectra of two sources observed at different times by the same detector. 
The same analysis procedure was used to obtain these inferred spectra. 
The Crab Nebula is a featureless relatively flat power law while 
Nova Monocerotis 1975 exhibits an equally featureless but steep thermal 
spectrum. 

IV. Discussion 

Cas A has no observable pulsed emission (Holt ^ 1973), and 

Fabian £t^ (1973) did not find a central intensity enhancement. Thus 

there is no evidence for a pulsar in the remnant. In the Crab Nebula 
the emission Is dominated by processes originating at the pulsar. This 
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featureless power law (Fig. 3) due to synchrotron radiation masks any 

possible thermal emi.ss^on from the remnant (but see Toor £t al. 1975). 

Without a strong pulsar, a complex thermal spectrum Is evident from 

Cas A. Charles £t (1975) fit their observation of Cas A to a two 

temperature thermal spectrum with kT * 2.6 keV and 0.7 keV. A thermal 

spectrum with temperature 1.4 keV was measured by Hill £t (1975) 

with increased silicon and decreased iron abundance. Serlemitsos ^ al . 

(1973) observed Cas A with a xenon detector similar to the OSO-8 CXS. 

They interpreted the 2-20 keV spectrum as a power law in photon number 

with a doppler broadened line feature centered around 6.9 keV due to 

+25 +26 

charge exchange between Fe and Fe cosmic rays and ambient inter- 
stellar matter. Their next most acceptable fit was a spectrum consisting 

of two thermal components. All report absorption consistent with 

22 -2 
N R: 10 atoms cm 
n 

In the present measurement improved statistics allow us to resolve 
this spectrum further. The feature around 6.7 keV persists, but is now 
seen to be narrow, broadened only to an extent consistent with the 
instrumental resolution at this energy (1.1 keV FWHM) . This can be 
interpreted as thermal emission lines from iron. Although the high 
energy (E > 7 keV) data is consistent with a steep power law (number 
index of 4.5), the presence of narrow iron lines strongly suggests a 
thermal source of temperature kT « 3.9 keV which also yields an accept- 
able fit for high energies. Around this temperature iron is most 
efficient in emitting K X-rays. For temperatures much below this, the 
highly ionized states are not populated and for higher temperatures 
the nuclei are normally fully stripped. However, in the 'standard* 
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supernova shock wave model presented in Gorensteln ^ (1974) (and 

references cited therein) the electrons and Ions In the plasma behind 
the shock wave may be at different temperatures. The time needed for 
the Ions and electrons to reach thermal equilibrium via coulomb collisions 
is long compared with the lifetime of the Cas A remnant. The electron 
temperature Is observed directly from the X-ray spectrum while the ion 
temperature is deduced from the velocity of the optical filaments and 
found to be ~ 30 times higher-about 40 keV. In thermal equilibrium at 
this temperc ture iron nuclei are -tripped of atomic electrons but in a 
non-steady state situation (cf. Kafatos and Tucker 1972) the ionization 
state could be less, leading to anomalous line emission. Alternatively, 
McKee (1974) proposed electrodynamic coupling as a mechanism to bring 
the ions and electrons into thermal equilibrium much tn^re i.uickly. rtt the 
electron temperature. Iron would emit Intense line emission In the 
steady state. In what follows we assume this last condition. 

The total energy flux In the iron lines can be compared to the 
continuum flux to obtain an estimate for Iron abundance In the plasma 
at this temperature. Several model calculations* have been performed 
for optically thin plasmas in this temperature range. Among them are 
Tucker 3»’.d Koren (TK) (1971), Felten, Rees, and Adams (FRA) (1972), 
and Raymond and Smith (RS) (1975). This last work Includes Gaunt 
factor calculations of Mewe (1972), line emission contributions from 
dlelectronlc recombinations, and new calculations of Ionization 
equilibrium. Table 1 displays results using the preceedlng models. 

The first two columns are inferred equivalent widths and detector- 
dependent line to continuum ratios respectively, with ’cosmic’ abundances 
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(l.e. logarithmic Iron abundance of 7.6 relative to hydrogen at 12.0.) 
Tlie third column shows the Iron abundance enhancement factor needed to 
produce the observed line to continuum ratio. 

The general conclusion we reach from Table 1 is that normal cosmic 
Iron abundance is adequate to explain this observation. If the lower 
ratio estimates from the first two models are correct, iron may be 
enhanced !n Cas A over cosmic abundances as Is to be expected If 
heavy elements are synthesized in highly evolved stars. In the early 
stage of the standard model the mass ejected by the supernova, 
is large compared to the Interstellar matter swept up, 

M ^ A n r^ m N (2) 

c 2 P ® 

assuming a spherical expanding shell into a homogeneous Interstellar 
medium of number density Nq. The radiating plasma matter will be 
more characteristic of evolved supernova constituents than of 
Interstellar matter. The shell luminosity is given by Gorenstein et al. 
(1974) as 

« 5x10^^ p(t) r 3 erg/s (3) 

2 

where P(T) is the emissivity of the plasma and R is measured 
in parsecs. In the present observation the 3.9 keV component has an 

>1 c _ 1 

intrinsic luminosity of 7.3x10 erg s and together with the 
known values of P(T) and R, we determine ■ 4.5 cm”^. This gives from 
equation (2),M^ i 4.6 M0. The value for the density is somewhat lower 
than previous estimates (Gorenstein et^ al. 1974) because only the flux 
fr»>n the high temperature component is considered. 
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If the higher estitnate for line to continuum ratio of KS is correct 
than iron is not enhanced in Cas A. Either iron was not over-abundant 
n the evolved star or Cas A is in a later stage than assumed in the 
preceeding and is less than the ambient mass swept up. However this 
poss^.bility does indicate that even cosmic abundance of iron should be 
readily observable by our and similar detectors via its line emission 
signature in hot optically thin gases. 

The low temperature component of the Cas A spectrum was 
observed as an excess in low energy channels over the high 

22 

temperature component. The observed kT . 7 keV and w 10 
Is similar to that in previous observations. Because of the 
uncertainties in the continuum contributions of the two components 
and the detector efficiency cutoff near 2 keV we are unable to 
estimate the silicon line contribution or comment on the need for 
enhanced silicon abundance. 

Supernova models with more than one thermal component have been 
applied to Cas A. McKee (1974) associates a high temperature thermal 
with the outward propagating supernova shock and a low temperature 
thermal with an inward propagating reverse shock. Recent results 
of Charles (1975) Indicate a spatial separation between components, 
with a low energy excess concentrated toward the center of the remnant. 
Chevalier (1975) attributes different components to hot gas as well 
as to the denser floculi moving away from the explosion site. Both 
models suggest a component with kT 20 keV. If this high temperature 
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componenc is present In addition to the observed 3.9 keV plasma, then 
between 3-50 keV the power in the former is at least a factor of 5 
below the power in the latter. 

We would like to thank the people involved in the GSFC experiment 
as well as the OSO project. In the construction and testing of 
the experiment we thank C. Glasser, F. Bitsa, T. Kaminski, .1. Saba, 

S. Jones, and E. Mundy. We also thank the delta launch crew at Cape 
Canaveral, and Hughes Aircraft Cc. for the excellent orbiting platform 
from which these observations were carried out. 
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Figure Captions 


1(a). The points and heavy bars are the observed counts In 
each energy bln with their statistical errors. The 
histogram reuslts from folding the best fit model 
described in the text through the detector response 
function. 

(b) . The inferred Incident spectrum obtained by dividing the 
data by the true channel-by-channel efficiencies. 

2. Residuals as a function of energy, obtained by subtracting 
from the data the best fit continuum folded through the 
detector response function. 

3. The spectra of (a) the Crab and (b) Nova Monocerotis 1975 

obtained with the same detector which observed Cas A 

but at different times. The vertical scale is photons 
cm~^sec”^keV”^ with different normalizations for the 


two spectra. 





600 k RESIDUALS ABOVE 



ENERGY (keV) 




(a) 

OSO-8 

CRAB 

1 1/6/75 



(b) 

OSO-8 

NOVA MONOCEROTIS 1975 
10/1/75 



10 

ENERGY (keV) 


100 I 


10 

ENERGY (keV) 


100 


